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Using product integrals we review the unambiguous mathematical representation of Wil- 
son line and Wilson loop operators, including their behavior under gauge transformations 
and the non-abelian Stokes theorem. Interesting consistency conditions among Wilson 
lines are also presented. 

Wilson loops provide a convenient method of obtaining gauge invariant observ- 
ables from lattices to strings. On the other hand the product integral formalism was 
developed in connection with matrix valued differential equationsEI, having built in 
the feature of order of the matrices. Consequently these come quite handy when 
dealing with path ordered quantities. As the theory of product integrals is well 
founded, we can use the various results in the subject to study Wilson lines and 
loops. In this paper we review the mathematics of the formalism of product integrals 
and its applications to the physics of Wilson lines and Wilson loops. 

1. Review of product integration 

The long history of product integration starts with Volterra§. In fact the product 
integral is to the product what the ordinary integral is to the sum. For starting 
point consider an evolution equation of the type Y'(s) = A(s)Y(x) on the interval 
[a, b]. Given Y(a) we want to find Y(6). We can obtain an approximate value 

Y(6) « e ^»H s "- s "-i) . . . e A ^ Sl - S0 ^Y{a) = H P (A)Y(a) (1) 

for any partition P = {so, s\, . . . , s n } of the interval [a, b]. In fact this approximation 
process can be made completely precise, and in the limit where the norm of the 
partition /i(-P) goes to zero, under suitable conditions, the expression (Q) converges 
to the product integral: Y(6) = linwp)_,.Q Hp(A)*Y(a). 

We briefly enumerate some of the basic properties of product integralsEIS. First 
we note that Yla e A ^ ds is the solution of the initial value problem we started 
with. From the formula det (J\* e M s ) ds ^ = e S a TlA ( s ) ds we m f er ^hat product 
integrals are non-singular matrices. We have a composition rule similar to that 
of ordinary integrals: J[ b a e A(s ^ ds = ]J b c e A ^ ds YF a e A ^ ds . The rule of differentia- 
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tion with respect to the endpoints reads: (jjy e A ^ ds J — A{x) Y[y e A ^ ds and 

The usual way of computing ordinary integrals is by the means of the Newton- 
Leibnitz formula, involving primitive functions. The analogous notion for product 
integrals uses the so-called L-operation. The L-derivative for a nonsingular matrix 
valued function P : [a,b] — » C„ xn is defined by: LP(x) :— P'(x)P^ 1 (x). It is 
immediate that Yl a e^ LP ' )l - s " >ds = P(b)P~ 1 (a), the analog of the well known Newton- 
Leibnitz formula. 

We reproduce three more important theorems that will be of great use later on: 
Sumrule: \[ x a e lM')+BW]de = p^) jjI e p ' 1 ^ B ^ p M ds ) where P(x) = Y\ x a e A ^ ds . 
Similarity rule: P(x) (HI e B ^ ds ) P'^a) = J]a e [iP(s)+P(s)B(s)P_1(s)lrfs - 
Derivative with respect to a parameter: For A(s; A) with the right propertiesa we 
have: £ ft x e A(s;X)ds = j* ds jj- e A( s -,x)dsdA( s . A) jj. eA(s , x)d ^ 

2. Wilson lines and Wilson loops 

Let M be an n-dimcnsional manifold representing the space-time and A a Lie- 
algebra valued connection one- form A(x) = A^(x) Tk dx^. For a curve C (closed) 

r b 

the Wilson line (loop) is defined as W\C\ = Ve*°- . If C is parametrized by a then 
A(x) — A((r)da, where A(a) = A^ (x(a))dx^ / 'da . A(a) plays the role of the matrix 
function in the product integral. 

Let us next consider the Wilson loop. For simplicity we assume that the loop may 
be taken to be the boundary of a two dimensional orientable surface £ parametrized 
by (a, t). We can construct the pulled-back field strength F a b in two different ways: 
as the the field strength of the pulled-back connection or as the pull-back of the 
field strength. It is easy to check that the two coincide. 

Our goal is to express the Wilson loop operator in terms of product integrals. As 
the Wilson loop depends only on the homotopy class of C, we can parameterize C 
conveniently: break it up into segments along which either a or r remains constant. 
The composition rule for product integrals ensures that the final result is indepen- 
dent of the intermediate points. Accordingly we have W[C] — W4 W3 W2 W\. For 
these intermediate Wilson lines cto = const, along W\ and W3, and a\ — const. 
along W% and W4. 

Based on the properties of product integrals, it was shown by two independent 
methodsEl that W[C] has a surface integral representation of the form: 

jx r TT P -y)F m (a' ;r')T(a' -y)da> ' dr' 

11 ° ' ( > 

TO 

where T is a given Wilson line. This can be interpreted as a non-abelian Stokes 
theorem. 

Let us now focus on the transformation properties of Wilson lines and Wilson 
loops under a gauge transformation A(x) — ► g(x)A(x)g~ 1 (x) — g(x)dg(x)~ l . It 
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was showna using product integration that Wilson lines have the expected transfor- 
mation Ve a). For the Wilson loop one has two ways to 

derive the transformation property: use the already established one for the lines, or 
use the surface inteeral representation provided by the non-abelian Stokes theorem 
(||). It was checkedn that both ways we get the well established answer, reinforcing 
our confidence in the techniques used. 

3. Further properties 

It is interesting to use the product integral formalism to derive interesting 
identities among Wilson lines. For brevity we introduce the following notation: 
P{A-b,a) = U b a e A{s)ds . Observing that f\ b a e A ^ ds = jfa e^M+t 1 -^ 8 )^ for 
any real A, based on the sum rule we can deduce immediately that 

b 

P(XA; a, b) P(A; 6, a) = J] e P(*A;a,s)(i-\)A( s )P(\A; S , a ) ds ^ (3) 

a 

An interesting immediate application of (^) is the addition law we obtain for A = ^ : 
P{2A; b, a) = P(A; b, a) Jjl e p ^ a ^ A ^ p( - A ^^ ds . 

Another application of (J3J) is by differentiating with respect to A. After some 
algebra, involving the formula for differentiation with respect to a parameter recalled 
earlier, we obtain the following result: 

j h a ds {P(XA; a, s)A(s)P(XA; s, b)P(A; b, a) - P(AA; s, b)P(A; b, s) 
[-P(XA; a, s)A(s) + (1 - A) J°dt [P(XA; a, s)A(s)P(XA; s, t)A(t)P(XA; t, s) 

+P{XA; a, t)A(t)P(XA; t, s)A(s)]P(A; s, a)]} = 0. (4) 

In conclusion we note that the results described here were extended to the s 
symmetric case, marking the beginning of supersymmetric product integratio 
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